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We present examples demonstrating that quasi-equilibrium fluctuation-dissipation behavior at 
short time differences is not a generic feature of systems with slow non-equilibrium dynamics. We 
analyze in detail the non-equilibrium fluctuation-dissipation ratio X(t,t w ) associated with a defect- 
pair observable in the Glauber-Ising spin chain. It turns out that 1/1 throughout the short-time 
regime and in particular X(t w ,i w ) = 3/4 for t w — » oo. The analysis is extended to observables 
detecting defects at a finite distance from each other, where similar violations of quasi-equilibrium 
behaviour are found. We discuss our results in the context of metastable states, which suggests 
that a violation of short-time quasi-equilibrium behavior could occur in general glassy systems for 
appropriately chosen observables. 

PACS numbers: 64.70.Pf, 05.20.-y; 05.70.Ln; 75.10.Hk 
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INTRODUCTION 

It is a common notion that in the non-equilibrium 
dynamics of glassy systems, fluctuations of microscopic 
quantities are essentially at equilibrium at short times, 
before the system displays aging. The concept of in- 
herent structures or metastable states 0], for instance, 
is closely related to this picture. There, one partitions 
phase space into basins of attraction, e.g., corresponding 
to energy minima. At short times the system is expected 
to remain trapped in such metastable states, leading to 
equilibrium-type fluctuations, but activated inter-basin 
transitions eventually produce a signature of the under- 
lying slow non-equilibrium evolution. 

A now standard procedure for characterizing how far 
a system is from equilibrium is provided by the non- 
equilibrium violation of the fluctuation-dissipation the- 
orem (FDT) 0j 0, • Starting from a non-equilibrium 
initial state prepared at time t — one considers, for 
t > t w > 0, 



OTw 



(1) 



where C(t,t w ) = (A(t)B(t w )) - (A(t)){B(t w )) is the con- 
nected two-time correlation function between some ob- 
servables A,B and R(t,t w ) = T 5{A(t))/5h(t w )\ h=0 is 
the conjugate response function; the perturbation as- 
sociated with the field h is 6TC = —h(t)B. Note that 
we have absorbed the temperature T into our definition 
of the response function. Thus, for equilibrium dynam- 
ics, FDT implies that the fluctuation-dissipation ratio 
(FDR) X(t,t w ) defined by Q always equals one. Cor- 
respondingly, a parametric fluctuation-dissipation (FD) 

plot of the susceptibility x(Mw) = L dr R(t,r) versus 
C(t,t) — C(t,t w ) has slope X = 1. In systems exhibit- 
ing slow dynamics and aging, on the other hand, the 
FDR X(t, t w ) is often found to have a non-trivial scaling 
form in the limit of long times [4J. Over short observa- 
tion intervals At = t — t w and at large waiting times i w , 
however, one expects to recover equilibrium FDT with 



X(t, t w ) — 1 according to, e.g., the inherent structure pic- 
ture. Such "quasi-equilibrium" behavior at short times 
has been observed in essentially all previous studies £| 
involving microscopic observables. 

Rigorous support is given to the quasi-equilibrium sce- 
nario in [fj. There a general class of systems governed 
by dissipative Langevin dynamics is considered. Based 
on the entropy production rate in non-equilibrium dy- 
namics, a bound is derived for the differential violation 
V{t,tJ) of FDT, 



V{t,U) = — 



R=[l-X(t,t w )} 



dC(t,t v 
dt w 



(2) 



The bound implies that under rather general assump- 
tions Q and for observables only depending on a finite 
number of degrees of freedom one has V(t, t w ) — > for 
any At > fixed and t w — > oo; we call this the "short- 
time regime" from now on. Via J3J) it is concluded that 
X(t, t w ) — > 1 in the short-time regime. Clearly, how- 
ever, the last step in this reasoning is only justified if 
\im tv/ ^ 00 (d/dt w )C(t,t w ) at fixed At is non- vanishing. 
This will be the case if, for example, C(t,t w ) admits a 
decomposition into stationary and aging parts Q at large 



C(t,t w ) — C s t(t — t w ) + C ag (t,t w ). 



(3) 



On the other hand, if C(t w , i w ) — ► for t w — > oo then also 
limt w _ >oo(9/c?£ w )(7(t, t w ) will normally vanish. For such 
observables, X(t,t w ) can a priori take arbitrary values 
in the short-time regime without violating the bound on 
V(t,t w ) derived in [5|]. 

In @, the stationary contribution can be defined for- 
mally as C st (At) = lim^-oo C(At + t w ,t w ). If C 8t (Ai) 
has a nonzero limit for At — > oo, then this is conven- 
tionally included in C ag instead. One then finds that 
the remainder, C ag , is an "aging" function: for large 
times it typically depends only on the ratio t/t w . Cor- 
relation functions of the form Q are often found in ag- 
ing systems, for instance when considering local spin- 
observables A = B = o~i in critical ferromagnets 0; here 
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C ag also contains an overall ^-dependent amplitude fac- 
tor. Further examples would be spin observables in p- 
spin models or density fluctuations in MCT [j|. In 
either case the stationary correlations C s t(t — t w ) are in- 
trinsically equilibrium-related [bulk fluctuations, dynam- 
ics within metastable states, "cage rattling"] and quasi- 
equilibrium behaviour is enforced by the bound of 
Beyond that, however, there are various examples in the 
recent literature [9|, llfl lll| displaying quasi-equilibrium 
behaviour even though C(i w , i w ) — * for t w — > oo. In the 
context of ferromagnets and spin-facilitated models such 
correlations are obtained when considering domain-wall 
or defect observables A = B = n*. A generalization of 
© accounting for the decrease of equal time correlations 
with i w is 

C(t,t w )=A(t w )c tt (t-t w ) + C B g(t,t w ). (4) 

In the short-time contribution C s t{t,t w ) = A(t w )c s t(t — 
i w ) H2, c s t is now defined as c s t(At) = lim^—xx, C(At + 
tw, t w )/A(t w ). It is useful to fix c st (0) = 1, so that A(i w ) 
must behave asymptotically as the equal-time correla- 
tion C(t w , t w ). Results for defect observables in one and 
two-dimensional Ising models are compatible with this 
scaling. For the one-dimensional FA and East models a 
similar picture emerges and it has been conjectured that 
aging contributions are in fact absent [13(101 • We discuss 
the scaling Q and its implications in more detail below. 
It is clear, however, that if C(t,t w ) ~ A(t w )c st (t — t w ) 
and also R(t,t w ) ~ A(t w )r st (t — t w ) in the short-time 
regime [where ~ denotes asymptotic similarity for t w — > 
oo] then quasi-equilibrium behaviour requires |33f that 
r s t(t-t w ) = (d/dt w )c st (t-t w ). The results of [jjllilEil 
indeed support this link. 

Recently, the authors of ^3 have exploited the notion 
of quasi-equilibrium to define a nominal system temper- 
ature Idyn even for models where a thermodynamic bath 
temperature T does not a priori exist, e.g., because the 
dynamics does not obey detailed balance. Tdyn ^ s deter- 
mined from the short-time dynamics of correlations and 
responses, and the authors of [l2| argue that for systems 
coupled to a heat bath this definition should generically 
reduce to Tc\ vr , = T. They in fact attempt a proof of this 
statement 34], for spin models in the universality class 
of the two-dimensional Ising model. 

The dependence of X(t,t w ) on the pair of observables 
A, B used to probe non-equilibrium FDT, Eq. l(T]). in 
finite-dimensional systems is still an actively debated is- 
sue 0,[3,[i3j[lJj[lM[i3- This is of particular relevance as 
regards the possibility of characterizing the slow dynam- 
ics in glassy systems by a time-scale dependent effective 
temperature T e g — T/X(t,t w ) ||. Beyond mean-field 
models one does not expect X(t,t w ), or more precisely 
its long time scaling, to be completely robust against 
the choice of A, B. Instead, it has been suggested that 
there may only be a limited class of "neutral" observ- 
ables which allow a measurement of effective tempera- 
tures 0, El • It then seems plausible that also the no- 
tion of quasi-equilibrium in the short-time regime may 



not hold for all observables. This prompts us to re- 
visit the observable dependence of short-time fluctuation- 
dissipation relations. 

To be able to carry out explicit calculations, we study a 
simple coarsening system, the one-dimensional ferromag- 
netic Ising spin chain with Glauber dynamics. Coarsen- 
ing systems are, of course, different from glasses but they 
do exhibit aging behavior, easily interpretable because of 
its link to a growing length scale. This makes them useful 
"laboratory" systems for testing general ideas and con- 
cepts developed for systems with slow dynamics. The 
dynamical length scale in a coarsening system - which in 
our case is just the typical domain size - allows one to 
distinguish equilibrated modes from slowly evolving non- 
equilibrium modes. For spatially localized observables one 
thus naively expects quasi-equilibrium dynamics at short 
times, as soon as the domain size has become much larger 
than the length scale probed by the observable. But, as 
we will see in the following, this is not true in general: 
there are many local observables that do not obey the 
equilibrium FDT even in the short-time regime. Such 
nontrivial violations of quasi-equilibrium behaviour have 
to be distinguished from what is found for global observ- 
ables such as, for example, the magnetization |9j: these 
depend on an extensive number of degrees of freedom, 
hence the bound of Q does not apply and one would 
not expect quasi-equilibrium behaviour at short times. 
In coarsening systems this is also physically transparent: 
global observables measure the dynamics on lengthscales 
larger than the typical domain size, where equilibration 
has not yet taken place. 

We analyze the non-equilibrium FDT in the Glauber- 
Ising chain for observables that probe correlations be- 
tween domain walls (defects) at distances d > 1. In Sec- 
tion [I] we define our two-time correlation and response 
functions; their exact derivation for the case d = 1 is 
sketched in the Appendix. We then recall some useful 
facts about the domain size distribution, both in and out- 
of equilibrium. In Sections ITT1 and IlIII we study adjacent 
defects, i.e. d = 1. Some features of the equilibrium dy- 
namics, where FDT is obviously satisfied, are discussed in 
Section[n] The low temperature coarsening dynamics are 
then analyzed in Section lTTTl and compared to the baseline 
provided by the equilibrium results. In particular, we fo- 
cus on the short-time regime in Sections 1111 AtlTlTCl The 
aging behaviour is discussed in Section IIII Dl while Sec- 
tion [^ni] deals with the crossover to equilibrium. Based 
on the understanding developed for d = 1, nonequilib- 
rium FD relations for defects at distances d > 1 are then 
studied in Section llVl We conclude in the final section 
with a summary and discussion. 



I. DEFECT PAIR OBSERVABLES 

In order to obtain nontrivial fluctuation-dissipation be- 
haviour in the short-time regime we have to consider 
non-standard observables; in the Glauber-Ising chain 
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local spin as well as defect observables satisfy quasi- 
equilibrium |qL Il7l Il8 j . However, as already mentioned 
in [9| , multi-defect observables are potentially interesting 
candidates for new results. The simplest choice in this 
class are the defect-pair observables A d = riirii+d with 
d > 1. We introduce the connected two-time autocorre- 
lation functions associated with Ad as 

C d (t,t w ) = (ni(t)ni + d(t)ni(t w )n i+d (t w )) 

-(n i {t)n i+d (t)){n i (t VJ )n i+ d{t v ,)). (5) 

The local two-time defect-pair response functions are 

5{ni{t)n i+ d{t)) 



Rd(t, t w ) — T 



Sh(t w ) 



(6) 



h=0 



where the perturbation STi. = —h(t)nini+ d is applied. 
Throughout the paper we use the short-hands C = C\ 
and R= Ri for the case d = 1. Our subsequent analysis 
of C and R is based on exact expressions, see Appendix. 
Since Ad is nonzero only if we simultaneously have defects 
at sites i and i + d, its behavior will reflect the domain 
size distribution in the system, and an understanding of 
the latter will be useful. 



A. Domain Size Distribution 

To summarize briefly, the Glauber-Ising chain is 
defined on a one-dimensional lattice of Ising spins <7j = 
±1 with Hamiltonian Ti = — JJ^i a i a i+ii where each 
spin <Ji flips with rate Wi(tr) = Ml — \~i(Ji{<Ji-\ + Cj+i)]; 
here 7 — tanh(2J/T). In terms of the domain- wall in- 
dicators or defect variables n, = ~ (1 — (Ji<j i+ i) € {0, 1} 
the density of domains Dk of given size k is expressed, 
using translational invariance, as 



Ac = (Ml - n x ) • •• (1 - n k -i)n k ). 



(J) 



As usual ( • ) refers to the ensemble average in the case 
of equilibrium and otherwise to an average over initial 
configurations and stochasticity in the dynamics. 

In equilibrium the derivation of D k for the Glauber- 
Ising chain is straightforward: from (Jli n i) — Ili^i) 
and translational invariance we have 



Dk, eq = {no} [1 - (no] 



ifc-i 



(8) 



The distribution of domain sizes k is thus exponential in 
equilibrium, with the most frequent domains those of size 
one. The mean domain size, on the other hand, is given 
by the inverse of the concentration of defects (no). One 
easily shows that 



Ceq(Teq) = (n ) = 



1 



cq 



V2j 



(9) 



For equilibrium quantities we generally use the equilibra- 
tion time T oq = 1/(1 — 7) ~ 5 exp(4 J/T) to parametrize 



temperature. At low T the defect concentration scales 
as c cq ~ (2T cq ) -1 / 2 . Hence the mean domain size is 
^(v^eq)- F rom ©: © the density of small domains 
with size k = 0(1) is then flat, D kt eq ~ l/(2r eq ). 

We note briefly that our Dk are densities of domains 
of given size, rather than a normalized domain size dis- 
tribution. The normalization factor is simply the defect 
concentration since, from (JSJ, J2T=i^ k ~ ( n o)- Ab- 
breviating c = (no), we can thus write Dk = cPk with 
Pk = 1- For small c, the normalized distribution 
Pk often assumes a scaling form, with k scaled by the 
mean domain size 1/c: Pk = cP(n) with k — kc, and 
correspondingly Dk = c 2 P(n). From (J8J, the equilibrium 
scaling function is exponential, P(n) = exp(— k). 

For the out-of-equilibrium case a derivation of Dk is 
rather less trivial; a corresponding calculation for the Id 
Potts model is given in |20|. For the Ising case and a 
quench at time t — from a random, uncorrelated initial 
state to zero temperature, the results are as follows: the 
mean domain size grows as c" 1 = O (yt) with typical 
domains having a concentration 0(c 2 ) = 0(t _1 ). For 
large domain sizes, k 3> yt, Dk(t) has an exponential tail 
O ((l/t)exp(-ak/Vi)); an expression for the constant a 
is given in |2(J • The density of small domains k <C \fi, on 
the other hand, is linearly related to the domain size, with 
Dk(t) — O (k/t 3 / 2 ). Correspondingly, the scaling form 
P(k) of the normalized domain size distribution decays 
exponentially for large k, but is linear in k for k <C 1. 

The precise scaling of the density of small domains 
k = 0(1) in nonequilibrium coarsening is easily derived. 
Instead of directly working out the Dk, which is cumber- 
some, consider for a moment the quantity 



D k = (n (l - 2rai) ■••(!- 2n k ~i)n k ). 



(10) 



In contrast to the Dk, any Dk can be conveniently ex- 
panded in terms of two-spin correlations. We have, in 
fact, D k = \{o x Ok - cfQOk - <xi<x fc+ i + (X a k+1 ). For zero- 
temperature coarsening one shows |'2l( 



D k (t) 



-21 



h(2t) 

2t 



(11) 



where the I n (x) are modified Bessel functions 22\. Now 
compare the definitions of Dk and Dk in the limit of 
large times. For both quantities we have that only 
states with uq = nk — 1 contribute. To leading order 
these states do not contain any further defects n, in the 
range i = 1, . . . , fc — 1, hence D k ~ D k ~ kf {Sy^rt 3 / 2 ) 
for t — > 00. States that do contain further defects in 
this range, on the other hand, cause Dk to differ from 
Dk- In an independent interval approximation, which 
gives the correct scaling but incorrect prefactors [2(], Hj| , 
the chances to have an additional defect at site i are 
DiDk-i = O (i~ 3 ) • Contributions from states contain- 
ing more than one defect in this range are even smaller, 
giving overall D k = D k + (i^ 3 ). By the same reasoning 
we also have (njTij+fc) = D k + O (t~ 3 )- These scalings 
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apply for any fixed k > 2 and in the limit of large t. 
For k = 1, finally, we have (njni+i) = Di = D\ as the 
definitions coincide in this case. 

Note that when comparing only the scale of typical 
domains in and out of equilibrium, an out-of-equilibrium 
system of age t is comparable to an equilibrium system 
with equilibration time r cq w t. Indeed, typical do- 
mains have size O (^/T eq ) and density O (t^ 1 ) in equi- 
librium, while out of equilibrium the same quantities 
scale like O {\f£) and 0(i _1 ), respectively. However, 
this correspondence does not extend to the details of the 
shape of the domain size distribution. In particular, it 
breaks down for small domains k — 0(1). In equilib- 
rium such domains have a concentration O (t^" 1 ) while 
in the corresponding coarsening situation their concen- 
tration O (fc/f 3//2 ) is much smaller. 

It is instructive to note that Glauber dynamics for the 
spin system \ corresponds to a diffusion limited reac- 
tion process [24J for the defects {rii}; the diffusion rate is 
|. At low T adjacent defects annihilate with rate close 
to one while pair creation, i.e., flipping a spin within a 
domain, occurs with rate l/(2r eq ). The latter process is 
important in equilibrium - continuously producing new 
domains of size one - but is unimportant at low temper- 
atures while the system is coarsening, and indeed strictly 
absent at zero temperature. This leads to the different 
scalings of the density of small domains in and out of 
equilibrium. 



II. EQUILIBRIUM 

In order to familiarize ourselves with the dynamics 
of defect pairs we now study the equilibrium behavior 
of the two-time correlation C(t,t w ). An exact expres- 
sion is obtained from the result for a quench to finite 
temperature given in the Appendix by taking the limit 
t w — ► oo at fixed At = t — i w . We use the notation 



C cq (Ai,T oq ) = lim tw ^ 
rium correlation; from 




) for the equilib- 



one has 



C eq (At, r cq ) 



1 



27Toq 
1 



H 1)Bq (At)|e- A *[/ -Ii](7At) 

--Hi,eq(At)l. (12) 



C<| 



Here we have introduced the short hand [-](x) to in- 
dicate that all functions enclosed in the square brack- 
ets have argument x. The function H is introduced in 
(IA.8I and discussed in the Appendix. Because FDT is 
satisfied in equilibrium the conjugate response to Ijl2|l 
is R eq (At, T cq ) = -dAtC eq ( At, T eq ) . Consequently the 
equilibrium susceptibility is given by 

Xeq(A<, r cq ) = C cq (0, T cq ) - C C q ( At , T eq ) , (13) 

and we subsequently focus on the discussion of C eq . 



A. Small At Regime 



Let us first consider the dynamics of C cq (At,T cq ) for 
finite At > and in the limit of low temperatures r cq 3> 
1. Via the definition JSJ of C = C\ the equal-time value 
C eq (0,T cq ) = (riin i+ i) - (nin i+ i) 2 is directly linked to 
the density of domains of size one, Di teq (r cq ) = 
From H, Pi and setting At = in (|T2|) 



-Dl,eq(T eq ) 



1 i-^/T^f 

2r fin T 2 



1 



-#l,eq(0). (14) 



At low temperatures Di teq (r eq ) and thus C eq (0,r eq ) 
scales as -Di je q( T eq) ~ l/(2r eq ). Now, for finite At > 
and in the limit of low temperatures r cq > 1 an expan- 
sion of l|12f> gives, to leading order, 

C cq (At, r cq ) = Pl (At)D lMq (T cq ) + O (r cq 2 ) , (15) 

where 



Pl (At) 



£ [/ 2 -J 2 ](At). 



(16) 



From (|15|l and our knowledge of the equilibrium do- 
main size distribution we may assign a direct physical 
meaning to pi(At): in the limit of low temperatures 
and at finite At the connected and disconnected cor- 
relations coincide to leading order, i.e., C eq (Ai,T cq ) ~ 
(nj(Af)nj+i(Af)rij(0)nj+i(0)). So only situations where 
sites i and i + 1 are occupied by defects at both times 
contribute to C cq . But since the size of typical domains 
scales as O (y/r^) , the probability for neighboring do- 
mains to be of size O(l) vanishes at low temperatures. 
Therefore, and since At = O(l), the defect pair at sites 
i, i + 1 at the later time At must in fact be the one that 
also occupied these sites at the reference time. Hence we 
may interpret p\(At) as the "random walk return prob- 
ability of an adjacent defect pair" . 

This scenario is easily verified by direct calculation. 
Consider a one-dimensional lattice containing exactly two 
defects at sites k and / with k < I at time t = 0. Denote 
by Pt(i,j) the probability to find these defects at sites 
i < j at time t. Since the T = dynamics of defects in 
the Glauber-Ising chain is diffusion-limited pair annihila- 
tion [2~ 

the pt(i,j) satisfy [2(| 



with diffusion rate ^ arid annihilation rate one 



d 1 

-Q t Pt{hj) = -2pt(i,j) + - \pt(i - 1, j) +Pt{i + 1,3) 



+Pt(i,j-1)+Pt(i,j + 1)] 



(17) 



This system of equations must be solved over i < j sub- 
ject to the boundary conditions pt(i,i) — 0. Using im- 
ages [2(1 El it is straightforward to show that the solu- 
tion is 



Pi 



(hj) = ^2Gt(i,j;k,l)p (k,l), 



(18) 



k<l 



where 



G t (i,j;M) = e~ 2t [h-klj-i ~ h-ilj-k] (t). (19) 
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It is clear from this result that the Green's function 
Gt(i,j; k, I) is in fact the conditional probability of find- 
ing the defect-pair at sites i < j at time t given that 
it was initially located at sites k < I. Consequently 
Pi(t) — Gt(i, i+1; i, as claimed above. The two-time 
defect-pair correlation C eq (At, r eq ), Eq. I|15|l. is thus to 
leading order given by the probability Di jeq (r eq ) of hav- 
ing a defect pair at sites i,i + 1 times the conditional 
probability pi(At) for this pair also to occupy the same 
sites a time At later. For At 3> 1 an expansion of 116|) 
gives Pi(At) ~ l/(27rA£ 2 ): the return probability for the 
defect pair drops quite rapidly as defects are likely to have 
disappeared via annihilation in the time interval [0, At] 
ifAi»l. 



B. Large Times 



lation C® q {At, r eq ) becomes Ai-independent. In other 
words, because of the rapid decay of the defect pair re- 
turn probability pi(Ai) we are more likely to find an 
independent defect pair at sites i, i + 1, rather than the 
original one, already on a time scale At — O (-/T eq ) . 
This is in marked contrast to spin or (single) defect ob- 
servables pi) , where this crossover happens on the time 
scale At = O (r cq ). 

Let us finally consider the equilibrium defect pair sus- 
ceptibility Xe q (A£,T cq ). According to it is strictly 
increasing, implying that R cq (At, r eq ) > at all times, 
and grows from its initial value of zero at At = 
to the asymptotic value C cq (0,r oq ) on an O(l) time 
scale. Explicitly we have from l|15[l the approximation 
X eq (A£,T cq ) L) l cq (r cq )[l - pi(At)] which holds uni- 
formly in At at low temperatures. 



When At becomes comparable to r eq the simple pic- 
ture discussed above breaks down; annihilation events 
with remote defects and pair creation are then relevant. 
But from equation (|12|l results for this regime, which are 
formally obtained by taking At, r cq — -> oo with their ra- 
tio fixed, are easily derived. In this limit we replace the 
modified Bessel functions appearing in (|12|) with their 
asymptotic expansions [2^|. This produces the leading 
order scalings 

C cq (At,r cq ) ~ 4n ^ t 2 Tcq [A<«r cq ], (20) 

C cq (At,r cq ) ~ J^L_e~ 2At /^ [At»r eq ]. (21) 

The expansion (|20|l matches the large At limit of Q15[l. 
So up to the time scale At = 0(r eq ) the decay of the 
connected two-time defect pair correlation is controlled 
by the defect pair return probability. For times At be- 
yond 0(r eq ), defect configurations are reshuffled via pair 
creation and the connected correlation vanishes exponen- 
tially as one might expect. For later reference we note 
that according to l[T5)l. (|2"U|l . (|2"T|) we have C eq (Ai, r cq ) > 
at all times. 

There is, however, a subtle effect in the underlying 
physics. This becomes obvious when considering dis- 
connected correlation functions. The disconnected de- 
fect pair correlation in equilibrium is C® q c (At,T cq ) — 
(ni(At)ni + i(At)nj(0)nj + i(0)), and is linked to the con- 
nected one via C° c (Ai,r cq ) = C cq (At,r cq ) + £> 2 , cq (T oq ). 
Now according to (fHjl we have -D 2 i0q (T oq ) = O (T e ~ 2 ) 
while from C cq (At,r cq ) = O (At^T^ 1 ) for At < 

T oq . Therefore, if At ^> y^^, C eq (At,T eq ) is negligible 
compared to D\ cq (T cq ) and so the disconnected corre- 



III. NON-EQUILIBRIUM 

In this section we analyze defect pair correlation and 
response functions for zero temperature coarsening dy- 
namics following a quench from a random, uncorrelated 
initial state. For the most part we will focus on the short- 
time behavior of these functions. Following our discus- 
sion in the introduction we decompose the two-time func- 
tions into short-time and aging contributions, 

C(t,t w ) = C st (Mw) + C ag (Mw), (22) 
R(t,t w ) — R st (t,t w ) + i? ag (t,t w ). (23) 

The two-time correlation and response functions are ob- 
tained from i|A.10|) and the construction of the response 
given in the Appendix. For mathematical simplicity we 
take the limit T — > but stress that the results are also 
valid for nonzero temperatures T > while the system 
is still far from equilibrium; as discussed in Sec. IIII El 
this requires t w ,t <gC r cq 2 / 3 . The response functions are 
always derived by taking the perturbing field h to zero 
before taking T — * 0. This has to be done to ensure lin- 
earity of the response: as discussed in more detail below, 
the size of the linear regime in the field strength h scales 
as T for low T. We find 

C 8t (Mw) = \e-^l%-lt]{t-tj-^, (24) 

flst(Mw) = \e- 2t [{ll-ll)-h{I -h)]{t-t w ) 
[I 1 + 2I 2 ](2t w ) 
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^ag (^) ^w) 



1 -2(t+t w ) ^1 ft + ^w) 

~4 (t + i w ) 2 



](2i w )- [la + h](t + t w ) 



h(2t w ) 
2^w 



1 -a(t +tw ) [Jq _ 2Ji + Ja](f _ t w )( [/l + ^ ] f + tw) [/ + /x](2M - [Jo + h](t + t w y 

■fie-^)^ - /,](* - t w) ( ^(^)/ 1 (H-t w ) _ [h-hKt + tw) [Io + /i](2tw) 



(26) 



Ji +2J 2 ](2t v 
It 



(27) 



r 



As we will see in Section IIII Al below, only the short- 
time functions Q24JI. (|25|l contain terms that contribute 
to leading order in the short-time regime. 

The results given above are exact. Before proceeding, 
we nevertheless compare them with simulation data to 
exclude the possibility of trivial algebraic errors in the 
derivation and confirm some of the more surprising fea- 
tures that are discussed below. For measuring two-time 
susceptibilities XA,B(t, f w ), with Aj , Bj generic local ob- 
servables, we use the standard method [25| of perturbing 
the system with STL = —h^2n eiBi for t > t w such that 
XA,B(t,t w ) = (T/h)J2 t h(Mt))} e in the limit h 0. 
Here the ti G { — 1,1} are independent, identically dis- 
tributed random variables and [-] e denotes an average 
over their distribution. Note that as our definition of 
the response contains a factor of T this is also the case 
for XA,B(t,t w ). Transition rates in the presence of the 
perturbation are, to linear order in h/T, 



.,(h) 



h 



(er) = w n (cr) + —w n (cr) [1 - w n (a)] ^ e «- 



x [B n - m (F n cr ) — B n _ m (<r ) 



(28) 



Here F n cr = (. . . , — a n , . . .) denotes the spin-flip opera- 
tor and w n (<r) are standard Glauber rates without per- 
turbation. In the procedure of measuring XA,B(t,t w ) 
the perturbing field h only appears in the combina- 
tion h = h/T. There is therefore a well defined zero- 
temperature limit at fixed h; the linear susceptibility is 
then obtained by using a sufficiently small h. We note 
finally that local perturbations only produce a few non- 
zero terms in the sum in (|28f) . The defect pair observ- 
ablcs Ai = Bi = nirii + i we are considering, for instance, 
only depend on <7j, Uj+i, and thus only m = 0, 1,2 
contribute to the sum in (|2"5|) . 

We show in Section llll Dl that with increasing i w it be- 
comes harder to see aging contributions in the two-time 
correlation and response functions. In order to be able to 
resolve these aging contributions in a simulation we have 
chosen t w = 1 and At = 0.1 ... 30. The data presented in 
Fig. ^ are for zero-temperature simulations in a system 
of 10 6 spins, averaged over 10 4 runs. For measuring the 
susceptibility we use h = 0.2, which is well within the lin- 
ear regime. Our code uses an event driven algorithm [26| 



although for the small times considered here a standard 
Monte-Carlo method would be just as efficient. As a full 
discussion of C(t,t w ) and R(t,t w ) will be given in the 
subsequent sections we comment only briefly on the data 
in Fig. ^ According to our decomposition (|2~2"|) we may 
think of the connected two-time defect pair correlation 
in Fig. Q[a) as the sum of short-time and aging contri- 
butions. From (|24|l short-time contributions are always 
positive while the aging contributions l|26|l are in fact 
negative, see Fig. ^b). From Fig. ^a) the correlation 
C(t,t w ) is dominated by short-time contributions up to 
about At — 1. In the range At = 1 ... 13 short-time and 
aging contributions compete, leading to a fast drop in 
C(t,t w ). At At ps 13, C(t,t w ) crosses zero, giving the 
cusp in the plot. For At > 13 the two-time correlation is 
negative, with short-time and aging contributions almost 
cancelling each other. The two-time defect pair suscep- 
tibility x(^w) shown in Fig. ^c), may similarly be re- 
garded as containing short-time and aging contributions 
according to (|2~3"|> . From Fig.^d) aging contributions in 
x(t, t w ) are tiny such that a plot of Xst(t, t w ) alone would 
fit the data shown in Fig. 2Jc) rather well. The small 
deviations of x(Mw) from Xst(Mw), that is Xag(Mw), 
however, are precisely as predicted by l|27(l . see Fig.QJd). 
Altogether, the data presented in Fig. are fully con- 
sistent with, and thus confirm, our exact results H24I27JI . 
We therefore now turn to a discussion of C(t, i w ), x{t, i w ) 
based directly on I|24I27|I . 



A. Short-Time Regime 

Here we analyze the dynamics of C(t,t w ) and R(t,t w ) 
in the short-time regime of At > fixed and i w — * oo. 
For the correlation we have from an expansion of i|26|) 
that the aging contribution scales as C ag (£, t w ) = O (t^ 3 ) 
in this limit; already at t w = 10 a plot of C ag would lie 
below the vertical range of FigQIb). The term C s t(t, i w ), 
Eq. on the other hand, is simply C st {t,t w ) = 

Pi(t — t w )Di(t w ) as follows from D\ = D\ and | (TT) |. 
IJTBJ. Because Di(t w ) = O (£ w ~ 3/2 ), aging contributions 
in C(t,t w ) are subdominant. In the short-time regime 
the connected two-time defect pair correlation is thus to 
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Figure 1: Comparison of our exact predictions IJ22H27II . shown as full curves, and simulation data [circles] at T = and i w = 1. 
Predictions for the susceptibilities are obtained by numerical integration of l|25|l and Q27ft. The data for aging contributions in 
panel (b) and (d) is obtained from the simulation data of C(i, t w ) and xtA^w) shown in (a), (c) by subtracting off the exact 
short-time contributions l|24|l [dashed line in panel (a)] and the susceptibility corresponding to l|25|l . In panel (b) and (c) error 
bars for the simulation data are invisible on the scale of the plot. The same is true for (a) and At < 10. For At > 10, on the 
other hand, the standard deviation of the data is around 10 -7 , which is reflected in the plot. In panel (d), finally, error bars 
[± standard deviation] are shown for every fifth data point but are no larger than the symbols. 



leading order given by C s t(t, i w ) alone, 

C(t, t w ) = Pl (t - t w )Dx(t w ) + O (t w 3 ) . (29) 

This scaling property is the key that separates short- 
time contributions C s t from aging terms C ag in (|22[l . 
Pi)l. Furthermore C st is of the form C at (t,t w ) = 

^4(*w)c s t(£ — t w ) as proposed in (0J. Comparison of fT5|> 
and H29|) shows that the nonequilibrium relaxation func- 
tion c s t (At) = Pi (At) coincides with its low-temperature 
equilibrium counterpart. Only the amplitude A is given 
by the dynamical density of defect pairs Di(t w ) instead 
of its equilibrium analogue Di iCq (T cq ). 

The result <|29[l is easily explained by random walk ar- 
guments, in full analogy to Sec. Ill Al Connected and 
disconnected correlations coincide to O (Df) — O (t~ 3 ) . 
For the disconnected correlation to be nonzero we need 
states containing a defect pair at sites i,i + 1 at time 
t w . These occur with probability D\(t w ) and to or- 
der O (t~ 3 ) there are no further defects in any finite 



neighbourhood. Hence only if these defects also occupy 
the same sites at time t is there a contribution to the 
disconnected correlation. This occurs with probability 
Pi(At) = G A t(0, 1; 0, 1) and ® follows. 

Now we turn to the behaviour of the defect pair re- 
sponse in the short-time regime. Expanding l|27|) shows 
that i? ag = 0(t- 3 ). But R st = 0(i„~ 3/2 ) from J2SJ) 
so R is dominated by i? st in the short-time regime, i.e., 
i?(t, t w ) = Rst(t,t w ) + O (iw 3 ) J us t as was the case for 
the correlations. We rearrange this result into the form 



R(t, i w ) 



ot w 



t-U 



xDi(t w ) + [t 



-5/2 



(30) 



where the expression in the square bracket coincides with 
the Ai-dependent factor in l|25|l . The ^-dependent am- 
plitude factor in l|25|l was expressed in terms of D\ — 
Di, Eq. C3J, using e- x I 2 {x) = e^I^x) + 0(x- 3 / 2 ). 
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Writing i? st (t,t w ) = A(t w )r st (t — t w ) as the equiva- 
lent of (J2J then shows that the noncquilibrium func- 
tion r st (At) is different from its equilibrium counter part 
(-d/dAt)c st (At) = (-d/dAt) Pl (At). In fact, from 
(d/dt w )C{t,t w ) ~ [(d/dt w ) Pl (t - t w )]Z>i(t w ) and (39 
we obtain for the FDR in the short-time regime, abbre- 
viating X{At) = lim^.^oo X(t, t w ), 



X(At) = - 



-2At 



(-d/dAt) Pl (At) 



At 



(31) 



Thus the FDR is neither equal to one nor even constant 
in the short-time regime. In particular, for At — > one 
finds 



X(0) 



lim X(t w ,t w ) = - 

t w —>OG 4 



(32) 



B. The Response Function 

Let us now try to understand the origin for the 
anomalous short-time response l|3U[l ■ For correlations we 
saw that the asymptotic equality between G and G st , 
Eq. (|29[) . could be easily explained by random walk ar- 
guments. This is, in fact, also possible for response func- 
tions in the short regime. We use that any response func- 
tion Ra,b as defined below Q can be written - via, e.g., 
the approach of [2l| - in the form 

RA, B {t,tw) = ^5ZA(0-')[pAt(0-'|F fe 0-)-pAt(0-V)] 
k er ,er' 

xA k B(a)w k (<r) [1 - w k {a-)]p u {<r). (33) 

This equation applies for general systems governed by 
heat-bath dynamics with Glauber rates Wk and for 
generic observables A, B. The p tw (er) denote state prob- 
abilities at time t w while PAt(""'|c) are conditional prob- 
abilities to go from state er to er' during the time interval 
At. In A fc S(<T) = B(F k er) - B{a) is just a short 

hand expressing the change of B under a spin-flip. 

In the concrete case of defect pair observables A = 
B = TiiTii+i we have AkB(cr) = except for k = i,i + 
l,i + 2. Denote the corresponding contributions to the 
response by 7*1 , r%, r% , respectively. Next work out a k — 
AkB{cr)wk{(T)[l — Wfc(er)]; it is convenient to use that 
Glauber rates at T = for the Ising chain may be written 
as Wk — + rife). It turns out that a 2 +i = and 

hence r 2 = 0. Also, r% and are related by reflection 
symmetry so we only discuss n. From a>i and (|33|l . 

r\ = - ^ fiim+i [pAt{er'\F l (T) - p At (a-'\a-)] 

x[rij_i(l - ni)n i+1 - (1 - ni_i)rijn. i+ i]p tw (cr). 

In the first line of this equation is to be read as rik(cr') 
while in the second one = rifc(cr). Now consider the 
term 7ii_i(l — rii)ni + i. It only contributes to r\ for states 
er containing defects at sites i — 1 and i + 1 but not 



at site i. For zero temperature coarsening and at large 
t w , however, we have that if there are defects on sites 
i — 1 and i + 1 then to 0(t~ 3 ) there will be no fur- 
ther defects in any finite neighbourhood anyway. We 
also know that the density of states containing defects at 
sites i - + 1 is {m-\rii + i) = D 2 (t w ) + O (t~ 3 ). Next, 
given any such state er, J2a-' n j(< T ') ri i+i(' T ')? , At(c r '|c) is 
the probability that these defects occupy sites i,i + 1 a 
time At later, that is Ga*(*, i + 1; i — 1, i + 1) from 119fl . 
The state i^er, on the other hand, has its defects on sites 
i, i + 1. [To see this, note that is a spin flip operator 
a.; -> -tTi, so n i _i(F i cr) = 1 - n i _ 1 (er) and n^Ficr) = 
1 - n t (cr) using n k = i(l - er fc cr fc+ i).] Consequently 
So-' n i( <T ') n i+i( <T ')PAt{ (T '\Bi (T ) = G At (i,i + + 1). 
Repeating this argument for the term (1 — ni_i)rijnj+i 
in n, where (n^i+i) = D 1 (t w ) + O (t~ 3 ), and working 
out r3 in the same fashion then produces 

R(t,t w ) - i[G At (0,l;0,l)-G At (0,l;0,2)] 

x [Di(t w ) + D 2 {t w )] + O (t~ 3 ) . (34) 

In this result invariance of G under index translation 
and Ga*(0, 1; — 1, 1) = G At (0, 1;0, 2) was used. From 
A(*w) = D 2{U) = D 2 (t w ) + O (t~ 3 ) and UU, 

(|19ll one shows that l|34|l precisely reproduces the short- 
time response R st , Eq. (|23|) . The subdominant O (t~ 3 ) 
corrections in l|34|) are aging contributions arising from 
multi-defect processes. 

The structure of l|34|l clearly reflects the mechanisms 
causing a short-time response. During the time inter- 
val [t w ,t w + 5t] where the perturbation STL — —hniUi+i 
is applied there is an increased likelihood for a defect 
pair located at sites i,i + 1 to stay there. The effect on 
(ni(t)n i+ i(t)) is accounted for by Ga*(0, 1; 0, l)Di(t w ) 
in (|34fl . Conversely, the chances for the defect on site 
i + 1 to move to i + 2 during the interval [t w , i w + St] are 
decreased. The corresponding change in (ni(t)ni + i(t)) 
is proportional to — Ga*(0, 1; 0, 2)Z?i(i w ). By the same 
reasoning and starting from configurations containing de- 
fects on sites i,i + 2 or i — 1, i + 1 one explains the re- 
maining terms in (|34l) . Overall, defects are on average 
closer to each other and more likely to occupy sites i,i + l 
at time t w + St due to the perturbation. However, this 
increases the chances for subsequent annihilation of the 
defect pair so that we should expect (m(t)m + i(t)) to 
become lower than without the perturbation eventually. 
Indeed, from ijMI) and D 2 (t w ) ~ 2Z?i(i w ) the instanta- 
neous response R(t w ,t w ) ~ |l?i(t w ) is positive. But 
as we increase At the response drops quickly and be- 
comes zero at At — t* m 2.132; here r* is the solution of 
G T - (0, 1; 0, 1) = G T * (0, 1; 0, 2). For At > r* the response 
is negative and ultimately vanishes as O (At~ 2 ) in the 
short-time regime. 

Our discussion so far explains the shape and origin of 
the short-time nonequilibrium response. But we still do 
not have an answer as to why R st differs from its equi- 
librium counterpart and thus violates quasi-equilibrium. 
To the contrary, from the above reasoning it actually 
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seems puzzling that we found i? eq (Ai,T oq ) > in equi- 
librium; see end of Section lllBI The answer to this prob- 
lem is non-trivial: although the rate for defect pair cre- 
ation is negligible at low temperatures, perturbations of 
this process contribute in leading order to the equilib- 
rium response. For coarsening dynamics, on the other 
hand, such processes are absent [at T = 0] or negligi- 
ble [at T > 0, see Section fill E| . Unfortunately, when 
using T > rates w n in l|33|) the simple random walk 
analysis from above cannot be repeated. We therefore 
limit our discussion to the instantaneous response. From 
PAt(c'\c) = da-' .ct at At = and setting A = B in 

Ra^w^w) = ^ i^kA( (T )] 2 w k (cr) [1 - w k (a-)]p K (cr). 

k,tr 

Writing Wk = |(1 — 7) + \^(jik-\ + Tik) and substitut- 
ing A = riiTii+i it is straightforward to show that the 
instantaneous defect pair response at T > is 



R(t w j tyy ) 



1 1 



■7 



2r t 



cq 



1 r 

2 



D 1 (t w )+ 1 2 D 2 {t w ) . (35) 



The first term in this result accounts for perturbations of 
the pair creation rate. Defect pair creation at sites i, i + 1 
corresponds to flipping spin c^+i within a domain, where 
Oj = = <7;+2- From the Ising Hamiltonian the asso- 
ciated cost in energy is A i+ iH(cr) = 4J. In the presence 
of the perturbation STi — —hniiii+i, however, this is low- 
ered to 4 J — h. Therefore, at T > 0, the perturbation 
increases the rate u>i+i of such spin flips and thus the 
density of defect pairs; recall that in the limit of low 
temperatures we first take h — > and then T — > so 
that the calculated response is always linear. From 1|35|) 
and at low temperatures, where 7 ~ 1, this produces a 
contribution of l/(2r eq ) in the instantaneous response. 
Now compare this to the other terms in (|35|) . using that 
D2 ~ D2 at low T whether in or out of equilibrium. 
For low temperature equilibrium we have -D2,cq ~ -Di.eq- 
Because -Di. oq ~ l/(2r eq ) we may write R eq ~ 2Di ieq , 
with all terms in l|35(l being of the same order. Out of 
equilibrium, on the other hand, D2 ~ 2D\. The first 
term in (|35|l is absent at T = or negligible compared 
to the others at T > for sufficiently small t w [see Sec- 
tion ^^Ej. Overall, we thus have the instantaneous re- 



sponse R r- 
equilibrium. 



D\ for coarsening but R c 



2D X , 



It is this difference in the prefactors that 



leads to X(t w , i w ) 



FD Limit Plot 



From (|29|l . (|30() we have that the two-time functions 
C(t,t-w), R(t,t-w) drop to an arbitrarily small fraction of 
their equal time values within the short-time regime, i.e. 
before they display aging. Therefore the exact FD-limit 
plot follows from the short-time expansions. Since the 
amplitudes of equal time quantities are time dependent 
we normalize x(Mw) = x(t,t w )/C(t,t) and C(t,t w ) — 




Figure 2: FD-plots for the defect pair observable Ai = 
mrii+i. The full curve shows the limit plot defined by 136L 
The dashed line has slope X = 4 and is tangent to the 
limit plot at the origin. The black dot on the right marks 
the end-point (l,Yi) of the limit plot. Simulation data are 
represented by circles. They are obtained from simulations 
at T = in a system of 10 6 spins. The probing field for the 
susceptibility is h — 0.5 and the data are averaged over 500 
runs. In the plot, t v — 10 is fixed and t = 10 ... 15 is the run- 
ning parameter. Inset: A plot of the FDR in the short-time 
regime X(At), Eq. 1311 . versus At. 



C(t,t w )/C(t,t) and plot % against 1 - C, see [S l27| . 
From J2SJ|, (1201) one obtains 



C = Pi(At), 

3 
4 



X 



l-pi(At)) - 2 / dre 



■At 



(36) 
(37) 



These equations apply in the limit f w — > 00 for arbitrary 
fixed At > 0. The resulting FD-plot is shown in Figure 
Note that when constructing FD-plots one generally has 
to keep t fixed and use i w as the curve parameter flq| . 
This convention ensures that the slope of the FD-plot 
is X(t,t w ). In the short-time regime we are exploring, 
however, the normalized functions only depend on t — t w 
and either t or t w may be used as the plot parameter. 
This is exact for t w — * 00 and correct to leading order in 
t w at finite t w . 

In Fig. |3 the slope of the plot at the origin, where 
t w = t, is given by X(0) = |. As At increases and reaches 
t* the response goes to zero. Consequently the suscep- 
tibility reaches a maximum at At — t* and the tan- 
gent to the FD-plot in Figure [21 becomes horizontal, with 
X(t*) — 0. As we increase At further the FDR (J3TJ turns 
negative and diverges linearly with At, X(At) — ► —00. 
Hence the FD-plot becomes vertical as it approaches its 
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end point (1, Yi). Taking At — > oo in (|o7[) . where the in- 
tegral is solvable , gives Fx = |- f « 0.545. Fluctua- 
tion dissipation relations for the aging case, where At and 
t w are comparable, are compressed into this point. So the 
plot in Figure |21 only reflects the fluctuation dissipation 
behaviour in the short-time regime. In order to demon- 
strate that the predicted violation of quasi-equilibrium 
can easily be observed in simulations we have included 
such data in Fig. 



D. Beyond Short Time Differences 

Our discussion of the non-equilibrium coarsening dy- 
namics so far was focused on the short-time regime where 
At > finite and t w — > oo; only the short-time terms 
C s t,i?st in our expressions (|22l) . 123|) for the connected 
two-time defect pair correlation and response functions 
contributed to leading order in t w . Let us now briefly 
summarize some interesting features of C(t,t w ) and 
R(t, t w ) beyond the short-time regime. Here At, t w — > oo 
simultaneously, and therefore the aging contributions 
C ag ,i?ag, Eqs. i|26|). must be taken into account. 

For correlations we expect to see an effect from the 
competition between the pair return probability p\(At) 
and the chance of finding an independent pair at sites 
i, i + 1 at the later time t, by analogy with the situation 
in equilibrium; see Section III Bl In non-equilibrium and 
for small At the disconnected and connected correlations 
coincide to leading order in t w . So from (|2"5|) the dis- 
connected correlation is C uc (t,t w ) « Di(t w )pi(t — t w ). 
Assuming that this equation applies up to sufficiently 
large At - though still much smaller than t w - we may 
now estimate the time-scale at which competition sets in. 
This is done by comparing C DC (£,i w ) to Di(t w )Di(t), 
which is the product of the independent probabilities 
of having a defect pair at sites i,i + 1 at time i w 
and at time t. Because we are assuming At <C t w , 
Di(f w )£>i(t) w D\(t w ). The scalings pi(Ai) = 0(At- 2 ) 
and £>i(i w ) = O (t w ~ 3/2 ) then show that C DC (t,t w ) be- 
comes comparable to Di(t w )Di(t) on the non-trivial time 
scale At = O (t w 3/4 ). In fact, from the plots in Fig. |3| 
the connected correlation becomes negative on that time 
scale. This means that for At t w 3 ' 4 , the chances of 
finding a defect pair at sites i, i + 1 at time t and at time 
t w are lower than those of independently finding pairs 
at both times: the presence of a defect pair at time t is 
negatively correlated with that at t w . 

We may picture this effect as follows. If we know that 
there is a defect pair at sites i, i + 1 at time t w , the neigh- 
boring defects are likely to be at a distance of the order 
of the typical domain size. Then, as time evolves, the 
original pair becomes more and more likely to have dis- 
appeared via annihilation while neighboring defects have 
not yet had enough time to reach sites i, i + 1. For the 
equilibrium domain size distribution these effects reach a 
balance on the time-scale At = 0(y/r^). In the coarsen- 
ing case, however, the relative concentration of small do- 




Figure 3: Top: Plots of the exact connected two-time de- 
fect pair correlation function C(t, i w ) for zero temperature 
coarsening dynamics, obtained from 12211 . 12411 . 1261 , The 
curves correspond to t w = 10 1 , 10 2 , . . . , 10 7 from top to bot- 
tom, respectively. The cusp in each curve separates the re- 
gion C(t, t w ) > [at small At] from C(t,t w ) < [at large 
At]. Bottom: Plots of the exact two-time defect pair suscep- 
tibility X;(t,t w ) obtained by numerical integration from 1231 . 

C?J. The curves correspond to t w = 10 1 , 10 2 , 10 3 , 10 4 
from top to bottom, respectively. 



mains - as compared to typical domains - is much lower 
than in equilibrium, so that annihilation of the original 
pair is comparatively the stronger effect. Thus, on the 
time-scale At = O (t w 3 ^ 4 ), we have a "hole" in the spa- 
tial distribution of defect pairs around sites + This 
hole persists up to the time-scale At = 0(i w ), where 
neighboring defects have had time to diffuse in eventually. 
The connected correlation function, see Fig. |3 therefore 
has three dynamical regimes at large t w : Up to times 
At <C i w 3 ^ 4 the expansion for the short-time regime l|29() 
applies and C(t,t w ) ~ C s t(t,t w ). In the time window 
t w 3 / 4 < Ai < t w the connected correlation is negative 
and ^independent, with contributions from C s t negligi- 
ble so that C(Mw) Df(t w ) (l/64)n-H- s . Fi- 
nally, at large At 3> t w the connected correlation remains 
negative but vanishes as C(t,t w ) ~ — (l/8)7r _1 Ai~ 3 , as 
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follows from expansions of (|24[1 . (|26[l . 

Comparison of the correlation functions in Fig. ^ and 
Fig. [3] shows that the simulation data at f w = 1 has 
only given us a glimpse on the full aging behaviour of 
C(t,i w ). From the scales of the plot in Fig. on the 
other hand, it is clear that such data are out of reach 
for simulations. Consider, e.g., the curve for i w = 10 3 : 
the connected correlation C(t,i w ) drops from its equal- 
time value C(t w ,i w ) ~ Di(t w ) of about 10 -6 to around 
10~ 12 , that is by six orders of magnitude, before it de- 
viates from its short-time behaviour C s t and displays 
aging. This illustrates the general problem associated 
with exploring the aging behaviour of correlation func- 
tions with a scaling of the form (0}. We have discussed 
this issue in the context of single defect observables in 
the Id and 2c? Ising models in There, long-time FD- 
plots are trivial with \ = 1 — C. However, this only 
reflects that quasi-equilibrium is satisfied and does not 
reveal any information about the aging regime. A sim- 
ilar situation is encountered in the lei FA model which, 
despite a trivial FD-plot [Tot fTl| . has X ^ 1 in the ag- 
ing regime j2^|. As regards the issue of measuring the 
asymptotic FDR X°° = lim tw ^oo lim^oo X(t, t w ), a so- 
lution for this problem was suggested in ;9j. It consists 
in using different observables which share the same X°° 
but are more easily accessed in simulations. 

The aging behaviour of the two-time defect pair re- 
sponse function R(t,t w ) is rather simple. Analysis of 
l(ZS). lf2"7|) shows that the short-time expansion R(t, i w ) ~ 
-Rst(Mw), Eq. lj3U|) . applies until At becomes compare- 
able to t w . More precisely, for 1 < < ( w , we 

have R(t,t w ) ( 3/32)7 r- 3 / 2 Af 2 i w ~ 3/2 . Note that 

as discussed in Sec. 1111 Bl the response is negative. In 
the opposite limit 1 < f w < At this crosses over to 
R(t,t w ) ~ -(9/32)7r- 3 / 2 At- 3 i w ~ 1 / 2 , accelerating the 
decrease of R(t,t w ) by a factor of t w /At. Intuitively 
speaking, the two defects that were located near sites 
i, i + 1 at time t w and caused the response are likely to 
have annihilated with other defects in the system when 
At 3> t w . This decreases the chances for such a defect 
pair to return to sites i, i + 1, and therefore the response. 

The scaling of R{t, t w ) has an interesting consequence 
for the susceptibility x(Mw) = f t drR(t,r), also re- 
ferred to as zero-field-cooled (ZFC) susceptibility. At 
large t this integral is dominated by the short-time re- 
sponse Rst, i-e. x{t,tw) ~ Xst{t,t w ). In the integral, as 
r departs from t the modulus of the integrand R(t, r) 
drops like (t — t)~ 2 for t — r ^> 1 and therefore the in- 
tegral converges quickly. Aging effects in the scaling of 
R{t,r) when t — r ^> r, i.e. r <C i/2, only give small 
corrections to the value of the integral. Therefore ag- 
ing contributions in the ZFC susceptibility x(^w) are 
subdominant. At large t and for any t w < t this implies 
X(i, tw) ~ Di(t)x(t — t w ) with x as given by i|57|l . Conse- 
quently we have x(^^w) ~ YiDi(t) in the aging regime, 
see Fig. [31 and contributions from aging effects in the 
response have vanished in x(^w)- We can see the ex- 
tent of this effect by looking back to Fig. ^ already at 



the very moderate value of t w — 1, aging contributions 
in x(^j^w) are marginal. While all of this is immaterial 
for exact calculations, where we start from the response 
in the first place, it is crucial for interpreting simulation 
results. In problems where the res pon se function has a 
scaling analogous to (@J , see 0, Hfl UJ$ , measurement of 
the ZFC susceptibility gives x(i,i w ) ~ Xst(t, t w ) and no 
information about the aging behaviour of R(t, t w ) can be 
extracted. In a measurement of the so-called thermore- 
manent (TRM) susceptibility p(t, t w ) = J Q W drR(t, r), on 
the other hand, this bias is not present since t— r > t—t w . 
So if At > i w , for example, the integral only contains con- 
tributions from R(t, t) with t — r > r and aging in R is 
revealed. The situation is precisely reversed as compared 
to the case of spin observables in critical coarsening [2!| , 
where the aging behaviour of R can be extracted from x 
but not from p. 

The nonequilibrium FDR X{t, t w ) as obtained from 
(|22|l - (|27f) has rather strange features when At and t w 
are simultaneously large. But since the observable Ai 
does not produce quasi-equilibrium FDT in the short- 
time regime we do not expect X(t,t w ) to have a sensi- 
ble meaning in the context of effective temperatures. We 
comment only that the short-time expansion (|31|l applies 
as long as At <C t w 3 ^ 4 , while the asymptotic FDR di- 
verges, X°° — oo. 



E. Equilibration 

Let us finally consider in more detail the crossover to 
equilibrium for a finite temperature quench T > 0. We 
focus again on the short-time regime and, for simplic- 
ity, discuss only the equal time FDR X(t w ,t w ). An ex- 
pression for (d/dt w )C(t, t w )| t=tw is obtained from (|A.10(1 . 
The instantaneous response R(t w , t w ) follows most conve- 
niently from (|35|l by working out Di(t w ) and Da(t w ) f Q1 ' 
a quench to T > [21| • The independence of the result- 
ing X(t w ,t w ) for three different temperatures is shown 
in Figure As expected the curves cross over from 
X(t w ,t w ) — | at sufficiently small i w [but still i w ^> 1] 
to the equilibrium value X(t w ,t w ) — 1 at large i w . The 
time scale for equilibration of A(i w ,i w ), however, is set 
by T cq 2 / 3 . In order to understand this result we con- 
sider the densities of small domains: in equilibrium we 
have the scaling -Di, C q(V C q) ~ l/(2T cq ) whereas for zero 
temperature coarsening D\{t w ) ~ 1/ (8 v^iw 3 ^ 2 ) ■ This 
shows that the dynamical density D\(t w ) is comparable 
to the equilibrium density -Di. cq (T cq ) for t w — O (r eq 2//3 ). 
By working out D\ = (riiUi+i) for a quench to T > 0, 
one easily verifies that D\ indeed becomes stationary at 
its equilibrium value for i w 3> T cq 2 / 3 . In this regime one 
also finds D2 ~ D\ as expected for low-T equilibrium, 
rather than D2 ~ 2Di in the coarsening regime. From 
the representation (|33|l of R(t w ,t w ) it then follows that 
the process of defect pair creation starts to contribute to 
the instantaneous response for t w » T cq 2 / 3 , and that the 
ratio of the other terms assumes its equilibrium value. 
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Figure 4: Plots of the equal time FDR X(t w ,£ w ) versus 
rescaled time t„/T ctl 2 ^ 3 for three different temperatures. The 
curves correspond to r cq = 10 3 (dotted), r cq = 10 4 (dashed) 
and T eq = 10 (dashed-dotted). The full curve applies in the 
scaling limit t w ,T cq — > oo with tw/Tcq 2//3 fixed. Inset: Plots 
of A^(t w , iw) versus tw for r oq = 10 3 ,10 4 ,10 5 . The hump at 
ty, = 1 is caused by transients following the quench. 



It is then not surprising that also the Ai-dependence of 
R(t,t w ) becomes identical to the equilibrium response 
throughout the short-time regime, i.e., quasi-equilibrium 



behavior is recovered for t w 3> r c 



2/3 



IV. NON-ADJACENT DEFECTS 

In the previous two sections we presented a compre- 
hensive discussion of the functions C(t,t w ) = Ci(t,t w ) 
and R(t,t w ) = Ri(t,t w ) for the observable A± — n,irii + i. 
Now we investigate to which extent our findings gener- 
alize to observables Ad = ni7ii +d which detect defects a 
distance d apart. As explained in the Appendix, an ex- 
act derivation of the associated functions Cd(t,t w ) and 
Rd{t,t w ) as defined in (J5J and (jfJJl, respectively, would be 
extremely cumbersome. Instead we exploit the fact that, 
in the short-time regime we are interested in, these func- 
tions can to leading order in t w be obtained from random 
walk arguments. 

Consider the states <x obtained from low temperature 
coarsening dynamics at large t w . In complete analogy to 
the case d = 1 we have that amongst the states er con- 
taining defects at sites i and i + d only a fraction O (t~ 3 ) 
have further defects in any finite neighbourhood. There- 
fore the discussion of (|15fl or l|29|) directly generalizes to 
any finite d > 1. In terms of the probability for a defect 
pair initially located at sites i,i + d to occupy the same 
sites a time At later, 



the correlaton C d (t,t w ) in the short-time regime is 

C d {t, t w ) = p d {At)D d {t w ) + O (t- 3 ) . (39) 

Although this argument applies directly to the discon- 
nected correlations it is also true for Cd(t,t w ) since both 
agree to O (-D^) = O (*w 3 ) m the short time regime. For 
response functions Rd{t,t w ) we use ()33fl and the same 
reasoning as in Section Till Bl to to obtain for d > 1 

Rd = + ~[G At (0,d;0,d)-G At (0,d;0 ) d-l)]A 1 _i 

+\ [Ga*(0, d; 0, d) ~ G At (0, d; 0, d - 1)] D d 

+\ [G A *(0, d; 0, d) ~ G At (0, d; 0, d + 1)] D d 

+ \ [G A *(0, d; 0, d) - G At (0, d; 0, d + 1)] D d+ i 



O ft 



(40) 



[In order to save space we have omitted here the time 
arguments Rd = Rd{t,t w ) and Dd = £*d(iw) etc.] An 
expression for the G's in 139I40|) is stated in l|19|) . giving 
in particular Pd(At) = e~ 2At [Iq — ij] (At), while we can 

estimate D d {t w ) = D d (t w ) + O (t~ 3 ) using JTTJl. 

Before we proceed with a discussion of |J35Jl, for 
nonequilibrium coarsening dynamics, let us briefly con- 
sider an equilibrium situation. The above assumption re- 
garding the nature of the states er that contribute to C d 
and Rd then still applies if the temperature is low. There- 
fore, to leading order in T cq , the equivalent of H15(l for 
d > 1 is C d , cq (At,T cq ) ~ pd(At)Dd, cq (T cq ) from J3§J|. In 
igDJl, on the other hand, we use D d± i )eq (r eq ) ~ D d , cq (T cq ) 
as the density of small domains is flat in low T equilib- 
rium. Combining terms then shows that Rd,c q {At, r oq ) ~ 
[(-d / dAt)p d {At)}Dd, cq {T cq ). Thus equilibrium FDT is 
recovered from (|39l) . (|40ll at low temperatures. This 
is non-trivial because we use zero temperature Glauber 
rates w n in the derivation of Q4Up. see Section Till Bl In 
contrast to the d = 1 case of adjacent defects, pair cre- 
ation processes do not contribute in leading order to the 
responses Rd with d > 1. This makes sense: the pertur- 
bation 5TL = —hriiTii+d acts on sites a distance d apart 
but pair creation is only possible on adjacent sites. So 
the pair creation rate for sites i + d, i + d+1 [say] is 
affected only if we already have a defect present at site 
i. The latter condition makes such contributions in Rd 
subdominant for d > 2. 

Having clarified this qualitative difference between the 
responses R\ and Rd with d > 2 we return to nonequi- 
librium coarsening dynamics. Here the density of small 
domains Dd(t w ) is to leading order proportional to the 
domain size d; more precisely we estimate Dd±i(t w ) us- 



ing JTTJ and e x I d±1 (x) = e x I d {x) 
allows us to rearrange gQj! into 



Ox 



- 3 /2). This 



Rd(t,t w ) 



t-t v 



Pd(Ai) = G At (i, % + d; i, i + d), 



(38) 



xD d (t w ) + 0(t- 5 ^ 



(41) 
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Clearly the nonuniform density of small domains pro- 
duces a nonequilibrium term in Rd- It therefore differs 
from its short-time equilibrium form. From l|39|) . (|41|) we 
obtain for the associated FDR in the short-time regime, 
again abbreviating X d (At) = lim^^oo X d (t, t w ), 



X d (At) = 1 - 



1 



-2At *d 



(-d/dAt) Pd (At) 



At ■ 



(42) 



The equations (|39[) . 1|41[) . 1|42|) have a surprising struc- 
tural similarity with their d = 1 counter-parts (|29|l , (f 301) , 
(|31|l . But, in contrast to the case of adjacent defect pairs, 
X d (0) = 1 for d > 2 rather than X x (0) = |, Eq.|23 From 
the discussion above we see that effects from perturb- 
ing the defect pair creation rate, which are subdominant 
when d > 2 but contribute to leading order for d = 1, 
are at the origin of this difference. However, it should be 
stressed that for At > all FDRs X d (At) deviate from 
unity on an O(l) time-scale. Therefore there is no quasi- 
equilibrium regime for any defect pair observable A d . In- 
stead we have from g2) that X d (At) = 1-0 (At 211 - 1 ) 
for At < d 2 while X d (At) = -Ai/(2d 2 ) + 0(1) for 
At > d 2 , see Fig.^l 

We visualize the violation of quasi-equilibrium be- 
haviour in terms of FD-plots. The scalings i|39|l and l|40|) 
show that with increasing t w plots become progressively 
dominated by the short time behaviour of C d and R d . 
A parameterization of the limit plots is obtained by tak- 
ing t w — ► oo at fixed At. Normalizing correlations and 
susceptibilities as in Section llll CI gives 



C d = Pd (At), 

Xd = (l-Pd(At)) 



A/ 



dre 



-2r I 2 d(r) 



(43) 
(44) 



The limit-plots for d — 2,3,4 are presented in Fig. 
Each plot has slope X d (0) — 1 at the origin [not shown] 
and follows the equilibrium line rather closely until At « 
0.25d 2 . Somewhere in the range 2d 2 < At < 2.5d 2 the 
plots reach a maximum where Xd = 0, and acquire a 
vertical tangent as they approach their end points (1,1^), 
with Xd{At) — > — oo diverging linearly with At. From 
lUD the Y d are 



Y d = l--(-l) d 
ird 



7T \ - 



fe=i 



(-i) fc 

2k- 1 



(45) 



giving Y 2 fa 0.962, F 3 ps 0.983, Y 4 ps 0.990 etc., and 
Yd = l — O (d~ 2 ) for large d. We remark that because the 
limit plots for d > 2 lie very close to the equilibrium line 
- Fig.[S]shows only the top-right corner of the plot - very 
accurate data would be needed to reproduce them in sim- 
ulations. Furthermore Xd(t,t w ) only converges slowly to 
its t w — > oo limit (|44|l . This is easily verified by numerical 
evaluation of Xd(t, t w ) = (1/C(t,t)) J drRd(t,T) based 
on 139|) . 140|) . These two facts in combination make it 
virtually impossible to see the limit plots, Fig. in sim- 
ulations. Equations (|39|) and l|4U|l . however, perfectly 




1-C, 



Figure 5: FD limit plots for defect pair observables Ad = 
rurii+d with d = 2, 3, 4 from bottom to top, respectively, ob- 
tained from 1431 1 . 1441 1. The dashed line has slope X = 1 
and represents equilibrium FDT. Black dots mark the end- 
points (l,Yd) of the limit plots. Inset: Plots of the FDRs 
in the short-time regime X d (At), Eq. 1421 . with d = 2, 3, 4 
from top to bottom, respectively, versus At/d 2 . The curves 
for d = 3, 4 are almost indistinguishable, showing that Xd 
quickly approaches a scaling form. 



reproduce simulation data for C d (t,t w ) and Xd{t,tw) 



J. drRd(t,T) already at small times, e.g. 
t = 10 ... 15 as used in Section lill CI 



10 and 



Our simple random walk analysis does not allow us to 
make predictions on the aging behaviour of Cd and Rd- If 
both At and t w are large, complicated multi-defect pro- 
cesses must be taken into account; only a calculation as 
sketched in the Appendix for the case d = 1 would allow 
one to study this regime. As regards the susceptibility 
Xd, however, we can predict that Xd{t,t w ) ~ YdDd{t) in 
the aging regime even though we do not know the pre- 
cise aging behaviour of Rd- This is simply because Xd 
is dominated by the short-time response as discussed in 
Section HTTT31 



Finally, for quenches to small but nonzero temperature 
quasi-equilibrium behaviour will be recovered when t w 3> 
T cq 2 / 3 , just as for adjacent defects. This follows from the 
dependence of the short-time response (|40|l on the density 
of small domains Dd and Dd±i and the fact that these 
densities level off towards their equilibrium values on that 
time-scale. 
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V. CONCLUSIONS 

In this paper we have analyzed the FDT behavior for 
defect-pair observables Ad = nirii + d in the Glauber-Ising 
chain. Contrary to the commonly held notion that short- 
time relaxation generally proceeds as if in equilibrium, 
none of these observables produce X = 1 in the short- 
time regime At <C t w ; this applies as long as t w is below 
the crossover timescale r oq 2 / 3 . We showed explicitly that 
this unusual behavior arises from the response functions, 
while the short-time decay of correlations does indeed 
have an equilibrium form apart from the expected over- 
all amplitude factor. The deviations of the responses 
from quasi-equilibrium behavior could be traced to two 
factors. First, in the out-of-equilibrium response for ad- 
jacent defects, events where pairs of domain walls are cre- 
ated are negligible, while in equilibrium they contribute 
at leading order. Second, all responses are sensitive to 
details of the domain-size distribution in the system, via 
their dependence on the density of small domains, and 
these details differ between the equilibrium and out-of- 
equilibrium situations. 

The inherent-structure picture mentioned in the intro- 
duction suggests a generic interpretation of our results: 
starting from an out-of-equilibrium configuration with a 
given number of defects or domain walls, we can loosely 
say that we remain within the same "basin" as long as 
no domain walls annihilate; the energy then remains con- 
stant. A similar interpretation has been advocated, in 
the context of the Fredrickson- Andersen model, in [3(| • 
Transitions to a basin with lower energy then correspond 
to the annihilation of two domain walls; at long coarsen- 
ing times f w , such transitions between basins are sepa- 
rated by long stretches of "intra-basin" motion. Within 
this picture, our defect-pair observables Ad measure pre- 
cisely when a transition to a new basin is about to hap- 
pen, i.e., they focus on the out-of-equilibrium, inter-basin 
dynamics. From this point of view it is not surprising 
that the Ad do not exhibit quasi-equilibrium behavior 
even at short times. Spin and single-defect observables, 
on the other hand, are not unusually sensitive to tran- 
sitions between basins, so that their short-time relax- 
ation is governed by the quasi-equilibrium, intra-basin 
motion [jj, fToL Hl| . This highlights the crucial depen- 
dence of nonequilibrium fluctuation-dissipation relations 
on the probing observables. 

The above interpretation suggests that lack of quasi- 
equilibrium behavior in the short-time regime could oc- 
cur quite generically in glassy systems. Certainly in glass 
models with kinetically constrained dynamics |3l| . one 
would expect observables that detect the proximity of 
two or more facilitating defects to display similar be- 
haviour to the one studied in this paper. More generally, 
the same should apply to observables which are sensitive 
to transitions between basins or metastable states. In 
structural glasses, conventional observables such as den- 
sity fluctuations are clearly not of this type. However, 
observables which measure, e.g., how close a local par- 



ticle configuration is to rearranging into a different lo- 
cal structure could be expected to display violations of 
quasi-equilibrium behavior. If the density of such con- 
figurations decreases with increasing t w then the bound 
of [j| is essentially void as discussed below Eq. @. It 
would be interesting to construct such observables ex- 
plicitly - thresholding of an appropriately defined free 
volume would seem an obvious candidate - and to test 
our hypothesis in simulations. 

Finally, the requirement that the short-time relax- 
ation should display quasi-equilibrium behavior could 
be used to narrow down the class of "neutral observ- 
ables" which are suitable for measuring a well-defined 
effective temperature in the limit of large time differ- 
ences. We note in this context that the condition 
limt w _» 00 (<9/<9i w )C'(£, t w ) ^ at fixed At is not necessary 
to obtain quasi-equilibrium FDT. For the single-defect 
observables considered in 0, IToL ITlj |. for instance, this 
limit vanishes yet quasi-equilibrium behavior is observed 
nevertheless. 
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Appendix 

We summarize below the ingredients that are needed 
to obtain our expressions for C = C\ and R = R\ based 
on the general results derived in The correlation 

C(t,t w ) is first reduced to multispin correlations by sub- 
stituting m = o (1 — fiCi+i) m ©■ This gives 

16C(f,t w ) = +C( , i),(o, i) + C(o,i), (1,2) — C(o,i), (0,2) 

+ C(1,2),(0,1) + C(l,2),(l,2) ~~ C(l,2),(0,2) 
— C(0,2),(0,l) — C(0,2),(l,2) + C(o,2),(0,2)i 

where 

-K(t)^(i))K(M%(*w)). (A.i) 

We have omitted the time-arguments in Cij in order 
to save space. When using the symmetries of Cij un- 
der translations, reflections and permutations [among the 
components of i and j but not between i and j] the above 
equation for C(t,t w ) assumes the simpler form 

C(t, t w ) = - [C^i^iUl) + C(l,2),(0,l) _ C(0,l),(0,2) 

_ C(0,2),(0,1)] + Y^C(0,2),(0,2)- (A. 2) 

Note that while C(t,t w ) can be expressed in terms of 4- 
spin correlations (|A.1I) only, this is not possible for d > 2; 
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in the latter case the expressions for Cd(t, t w ) contain 8- 
spin two-time correlations and exact calculations become 
exceedingly cumbersome. In full analogy to C(t, i w ) the 
response R may be decomposed into 



1 



-R(Mw) — o [-^(o,i), (0,1) + R 



(1,2), (0,1) 



R 



1 



-#(0,2), (0,1)] + :[g-R(0,2),(0,2)) 



with 



Ri,j — T 



5hj(t v 



(0,1), (0,2) 

(A.3) 
(A.4) 



For the multispin response function (|A.4J) the field hj is 
thermodynamically conjugate to Gj l Uj 2 . Let us remark 
that while the pair C, R violates quasi-equilibrium this 



';3> Ri,J- 



is not the case 35] for the constituting pairs Cij 
Next, the 4-spin correlations in (|A.2() are expressed in 
terms of the result given in |2l| . viz., 



+ [Km - H ^-n ( 2Ai > 2 M] a n -n (2t w ) 

-[ + e- At I il . h (lM) + £l h ] 
x[-e- At / l2 _ J2 ( 7 Ai)+£/ 2j2 ] 



-A/ 



-A/ 



/ lW2 ( 7 At)+£/; j2 ] 



(A.5) 



Here and below the indices i = (11,^2) and j = (jx,j2) 
must satisfy ii < 12 and ji < j2- The multispin response 
function I) A. 411 for the case j 2 = Jx + 1 is also stated 
explicitly in |21| and reads [j ,s = (ji — and j 2 ' s = 
U2J2 + I)] 



Ri 



-e~ At I 



-A/ 



-21 



-e- At 7, 



H-22 



-e _A *7 



r 
2 



± ) (- e - A *7 fe _,,+£ 



2 

2f 
2 

9 

2 



12 -32 ^ j 2 



-p~ A * 7 ■ + £ j 



-p _At 7 ■ + f j 



-e _At 7 ■ + £ j 



7_ 
2 



-At 



-At 



-fsi-Ji+1 + ^1,31-1 



-At/". . 

J l2-J2-1 T C l 2 J 2 + l 



e- A *7, 



J1-J2 



^2,a 

- 1 ~T °ilj2+l 



(A.6) 



We have omitted the arguments (7 At) of the functions 7 in order to save space. The multispin response functions for 
the case ji > j\ + 1 are not given explicitly in |2l| . However, by following the general procedure developed there one 
verifies the result \j l,s = {ji - 1, jx, ji + 1, j 2 ) and j 2 ' s = (J1J2 ~ 1, hi h + 1)] 



n _ 1 „-At r 
iti.j T c 1 i\— ji 



7 



— IX — i— 1 |-e- At 7 



-(+e- At 7 



12-31- 



+e- At 7, 



%2— 31 



2 - 



7 



l-^l(- e - A *7 



(+e- A %_ jW + ££; i+1 ) i^ 2 _ J1+1 + ( < v /, „ +e£ a ) II 



' £>2,i2 


)4[ 


■ + (" 






>-£[ 


■ + (- 


e- At 7 4l 



z j2-ii-i 



^ _ 77 



e At 7 £l _, 1+ i + £ 



»i-ii+i " | " c, »i,ji-i J ^j'a-Ji-i 



-e~ At 7 ■ 



{+ (1 - y) (+^ At + Ci) - y [ (+^ At 



(-e At 7i 2 _j 2+ i + £ i2 ,j2-i) H 32-]i+i + (~ e At ^2-ji-i + j2 +i) Z'.-.' i 



e- m /i,-,-^ - ( 1 - 4- I ( I n - n + ^ 



,.)4[(- 



■e- Ai 7 il _ il +^X)77 2 



^-e ^J^-jj+x + H j2 -j 1+ i + (-e A *7 il _ j2 _ 1 +^,j 2+1 ) -^2-31-1 j 



(A.7) 



Again, all functions 7 have argument (7At) and addi- ment (2t w ). After substituting (|A.5|I . (| A. 6|) . (| A. T|> for 
tionally all functions H appearing in (|A.7|I have argu- 
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the multispin correlation and response functions in IIA.2I) 
and ljA.31) . C(t, t w ) and R(t, i w ) are expressed in terms of 
/, H, £ and T . Next we also represent £ and T in terms of 
/ and H by applying the corresponding formulas derived 
in H3, viz. 

= H jly _ ie (At,2t w ) 

- S 3> m e ~ At he-mil M) fl J -„_ m (2t w ), 

and 

^.i. = H l5 _ u (2At,2t w ) 

- Sj,m e- A */ le _ m ( 7 At) H ie - m (At, 2t w ) 

m 

+ 5 i>™ e~ At Irs- m (lAt) H ie - m (At, 2t w ) 

m 

+ ^ 5 s , m S jtn e- 2At J is _ ro ( 7 A*) / i5 _„( 7 At) 

X -^n — m ( ) ; 

where 5j iTn = 1 — sgn(ji — m) ■ • ■ sgn(jp — m) with 
p = dim(j), i.e. p = 2 or 4 in our case. The sums in these 
equations are finite since 8j m is nonzero only within the 
index-range covered by the components of j. Upon sub- 
stitution of the latter equations for all functions £, T 
the defect-pair correlation and response functions are ex- 
pressed purely in terms of I, H. The functions H , in 



turn, are expressed in terms of modified Bessel functions 
via 

fl»(ti,f 2 ) = | J dTe~ T [/„_! - J n+ i] ( 7 r), (A.8) 

where we use the notation H n {r) = H n (0,T). In equi- 
librium the quantity 7J„ i0q (r) = limt_,oo H n (r, t) is rele- 
vant, cf. Eq. Simplification of the expressions for C 
and R are possible when using the recursion formula 

2 

H n +i{ti,t2) = — H n -\{t\, £2) H — H n (ti,t2) 

7 

+e- 41 -* 2 [J n _i-/ n+ i](7(ti+t 2 )) 
-e- (l [/„_! - / n+1 ] (Ttx) . (A.9) 

One easily proves (|A.9|I when substituting l|A.8(l and inte- 
grating by parts. With i?o(ii, ^2) = 0, which follows triv- 
ially from HA.8|I . any function H n {ti,t2) may thus be de- 
composed into modified Bessel functions and H\(t\,t2)- 
Also, the recursion \I n -\ — I n +i](x) = ^I n (x) is useful 
for rearranging the results. We use Mathematica 5.0 to 
carry out the algebraic manipulations described above. 
The procedure yields significant cancellations in the ex- 
pressions for C and R. For a quench to T > we obtain 
see below, and a similar expression for R; taking 
T — > where the integral Hi is soluble [2l| then produces 
the results (j2"2"|) - (|2"T)) . 



C(t,tw) 



1 



4 7 2 

+Hi(t-t w ,2t 
1 



(t+tw) J i(7^ + M) 
2 



(t+tw) /i(7(* + *w)) _ 27e -(t-t w ) [7o _ /i](7(f _ fw))[1 _ Hi{2tw)] 



t + t 



• eq 



(t+M /i(7(* + *w)) _ 27e - ( ^ w) [/o _ 7i] (7(t _ u)) fj_ + h^y\ 



— ^ w > ^ v 



(A.10) 
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